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Here we show that, in principle it is possible to clone (measure) a single arbitrary unknown 
quantum state of a spin-| particle (an electron) with arbitrary precision and with success probability 
tending to one, using protective measurement. We first transfer the information from spin to spatial 
degree of freedom (d.o.f) of system electron, then trap it in a double well potential, and finally 
measure it protectively using a probe electron (which donot get entangled with system electron, but 
still extracts expectation value of an observable from a single quantum system (system electron)) 
to obtain information about the unknown spin polarization. Nonorthogonal state discrimination 
being a subclass of cloning, part of the paper (till finding out 0^, polar angle corresponding to the 
unknown spin polarization) is sufficient for discrimination. 


I. MOTIVATION 

“ There is no law other than the law that there is no 
law...All is mutable J A Wheeler [1]. Then how about 
the ‘no-cloning’ law? Assume there exists a unitary op¬ 
erator U which can clone two non-orthogonal states as 
follows: 

C/|0)|0) = |0)|0), U\+)\0) = |+)|+) (1) 

where, |-|-) = (|0) -I- |l))/-\/2, and, {|0), |1)} are eigenkets 
of z-component of total spin angular momentum (Sz). 
Now taking the inner product of two equatons in 1, we get 
1 = y/2, which is wrong. Hence, generalizing, we can say, 
it is impossible to clone an arbitrary unknown quantum 
state (which is non-orthogonal to the reference state |0), 
in eq. 1, of course excluding |0)) via a single unitary 
operation, which is the essence of no-cloning theorem. 
Now instead consider the following operation: 

??|0)|0) = |0)|0), D|+)|0) = |+)|+), (2) 
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where, D is an arbitrary linear operator on the Hilbert 
space, 'H 2 0 ^2, of two spin-1/2 particles. Operators 
{|crfe (Si <71] k, I = 0,1,2, 3}, where, (Jq = 1 (2 x 2 iden¬ 
tity matrix) and ai’s {i yf 0) are Pauli matrices, form an 
orthonormal basis in the product Liouville space L ® L 
of operators on "^2 ® ^2 [2] . hu are the coefficients of 
decomposition of D in the basis { 5Cfe 0 trz}. Therefore, 
we have 16 unknown complex coefficients h^i i.e., 32 real 
coefficients and 16 constraint equations (8 real -I- 8 imag¬ 
inary, obtained from eq.s in 2). Hence we can express 16 
real coefficients in terms of 16 other arbitrary real param¬ 
eters. This shows there exists infinitely many solutions to 
eqations in 2 (for one such solution see [3]). D being sum 
of unitary operators, CTfc 0 cr;, is a non-unitary operator in 
general and specifically in eq. 2. It may not be possible 
to directly implement the operator D in an experiment 
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(however, there might be some indirect way). As there 
are infinitely many solutions to eq.s in 2, many of them 
may be just mathematical objects with no relevance to 
physical operations. Now let us generalize the case in eq. 
2 : 

i^|0)|0) = |0)|0), i^|m)|0) = |m)|m), (4) 
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where, Irh.) (given by eq. 8) is the unknown state to 
be cloned. Again there exists infinitely many solutions 
to eq.s in 4 (for one such solution see [4]). In eq. 4, 
if |m) = |1) then, T>|I)|0) = |1)|1). But, D is a linear 
operator. Hence, 

_D|m)|0) = cos ^??|0)|0) -I- sin ^e*'^’"Zl|l)|0) (6) 

= cos^|0)|0) -bsin^e*'^’"|I)|l) 7^ |m)|TO) (7) 

Hence, if D is linear, we cannot clone an arbitrary un¬ 
known state [5]. Hence, D must be a nonlinear operator 
(i.e., eq. 6 doesnot hold) apart from being nonunitary, to 
clone an arbitrary unknown quantum state. The myste¬ 
rious quantum measurement process, which involves am¬ 
plifying information to the classical limit, is one such 
operation which is both nonlinear [6] [7] and nonunitary 
(as information is irreversibly amplified to the classical 
limit which increases entropy, eg., electron absorbed by 
the screen in Stern-Gerlach experiment). Hence, it mo¬ 
tivates us to ask the following question: ‘Is it possible 
to clone an arbitrary unknown quantum state, through a 
combination of linear unitary and nonlinear nonunitary 
(measurement) operations ?’. We found the answer to be 
‘Yes’. In the protocol that we are going to describe, we 
carry out many such measurements (nonlinear nonuni¬ 
tary) eg., on probe electrons (which donot get entan¬ 
gled with system electron, but still extracts information 
from it!), there by relaxing both unitarity and linearity 
constraints, which allows us to clone. In the following 
protocol, even though operations on system electron are 
linear unitary, operations on probe electrons and other 
auxiliary systems are nonlinear nonunitary, and clubbing 
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informations obtained via measurements (which is some 
kind of sum of unitaries which is nonunitary in general) to 
prepare another qubit in the state same as that of given 
qubit, and to get a final state similar to that on RHS of 
eq. 4, makes the global (system+probes+observer) pro¬ 
cess nonlinear nonunitary. 

Finally we like to mention that, in the cloning protocol 
that we are going to describe, in principle we can clone 
a single arbitrary unknown quantum state of a spin- 1/2 
particle (electron) with success probability tending to one 
(but not exactly equal to one, and hence there is probabil¬ 
ity, however small, of failure) and with arbitrary precision 
(which, strictly speaking, does not correspond to an ex¬ 
act or identical copy, as there is an error, however small). 
Hence, strictly speaking, it is not perfect cloning (exactly 
zero error) with success probability exactly equal to one. 


II. CLONING PROTOCOL 

Protective Measurement was invented by Y Aharonov 
and L Vaidman in their seminal paper [ 8 ]. The cloning 
protocol that we are going to describe is motivated by an 
idea proposed by Y Aharonov, J Anandan and L Vaid¬ 
man in [9]. Let the system be an electron which is a 
spin-i particle with charge — e and mass M. Let its spin 
be in an arbitrary unknown state: 

|to) = cos^lO)-bsin^e*‘^™|l) (8) 

where m is a unit vector along unknown spin polarization 
given by, 

m = sin 6m cos (pmi + sin 9m sin (pmj + cos 9mk (9) 

where, i,j, k are unit vectors along positive x,y,z axes 
respectively of fixed lab frame, dm and (j)m are polar and 
azimuthal angles respectively. For convenience, let the 
system electron’s wave packet be a Gaussian centered 
at zo = 0 , poz(average momentum)= 0 , which satis¬ 
fies AzApz = h/2, say, at time t = 0. Az and Ap^ 
are standard deviations in position and momentum re¬ 
spectively. The variance Apl, can be Squeezed (at the 
cost of increasing Az^) and made small enough (but not 
zero, because, if zero then we have a plane wave, which is 
not square integrable and hence, rigorously, cannot rep¬ 
resent a physical state of system electron (p23 [10])) to 
be taken as zero for all practical purposes. For eg., if Az 
is of the order of 10 “^m then Ap^ will be of the order 
of 10~^^kg m s~^, which can be treated as zero in con¬ 
junction with eq. 29, for all practical purposes. Also, in 
principle we can make Apz arbitrarily close to zero (but 
not exactly zero, there by satisfying square integrabil- 
ity requirement). To compensate for the corresponding 
increase in Az, in principle we can make width ‘a’ of po¬ 
tential well (Fig.l), arbitrarily large as explained below 
eq. 31. Hence, in principle we can push the system elec¬ 
tron into one of the eigenstates of i?os( 12 ) or Hs (18), 


as discussed below eq.s 14 and 28 respectively. 

As the wave packet evolves, it spreads in position space 
but not in momentum space. Its average momentum 
(poz) as well as momentum dispersion Ap^, do not change 
with time. Hence, at t > 0, AzApz > fi/2. Because 
Pqz = 0, even zq (center of wave packet in position space) 
do not change with time (pp 64-65 of [10]). Now an in¬ 
homogeneous magnetic field along z-axis is switched on 
for an interval of time r. System electron evolves under 
the total Hamiltonian p^/(2M) -|- Hint, where, interac¬ 
tion Hamiltonian, Hint = —p ■ Bizk = —'jBiSzZ, Bi is 
the constant gradient(Tes/a m~^) in magnetic field along 
z-axis [11], 7 is the gyromagnetic ratio of electron which 
is negative. As poz = 0, in the limit Apz —>■ 0 (case of 
our interest, see below eq. 31), p1/{2M) —>• 0. However, 
as it evolves under Hint, it gains kinetic energy. We can 
choose static field Bq [11] such that actual interaction 
energy H[^^ » p^/(2M) for any Bi, r [12]. Hence, we 
can neglect the evolution under p^/(2M), and write: 

IV’(aPz)) = exp(--i?™tT)lm)x(pz - 0) = 

cos^|0)x(pz -7 S*^t) -ksin^e*‘^™ll)x(pz -byS^^r) 

( 10 ) 

where x{Pz ~ 0) is the wave packet of system electron 
in momentum space centered at poz = 0 (see eq. 40 for 
definition of wave packet), which has split and got entan¬ 
gled with its spin d.o.f. Taking its inverse Fourier trans¬ 
form and finding its time evolution, we see that they are 
still entangled and the split packets are moving in oppo¬ 
site directions. Gradient Bi can be made either negative 
(p388 of [10]) or positive. Let Bi < 0. Then, system 
electron is in a superposition state of having momentum 
'jBi^T (hence moving along positive z-axis) with prob¬ 
ability cos^ ^ and having momentum —^Bi^r (hence 
moving along negative z-axis) with probability sin^ 
Let’s trap the system electron in a symmetrical double 
well potential (trapping procedure explained below eq. 
31) with the following description(see for eg. ,[10, 13]), 

f 0, for (6 — a/2) < z < {b + a/2) 

Vi{z) = < 0, for -{b + a/2) < z <-{b - a/2) (11) 

I oo, everywhere else 

where, a is the length of each potential well and b is the 
distance from origin of coordinate system to center of po¬ 
tential well(see Fig.l). Vi{z) is the potential experienced 
by system electron only, but not by any other charges 
outside the potential well. 
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FIG. 1. Symmetrical infinite donble well potential 


and hence its average kinetic energy is ^ = 
, where all the parameters are known. We 
can choose the values of free parameters a, 6, Bi and r (as 
explained below eq. 14), such that the system electron 
enters one of the eigenstates of the system Hamiltonian: 

H.s = i^ + Vr{z) ( 12 ) 

where, Vi{z) is as given in eq. 11. State of the system 
electron after trapping it in double well potential is given 

by, 


System electron has gained an average momentum 


l^ 0 sn( 2 :)) =cos^|0) e{z-{b-^))-e{z-{b+^)) y^sin(A:o^„(6 +“-z)) 


Or, 


■ Sin 


e*^-|l) B{z+{b+'^))-9{z + ib-^)) 


^ sin {kosnib+ - + z)) 


(13) 


where, fcosn = n = 1,2,3,.... Unit step function is 
defined as. 


9{z - zo) 


0 , for z < Zo 
1 , for z > Zo 


(14) 


which implies ^6*(z — zq) = S{z — zq), Dirac-Delta func¬ 
tion. We note that, even though Hqs commutes with 
parity operator, its eigenstate |4’osn(2;)) is not symmetric 
due to entanglement with spin d.o.f. Similarly |'l>oara(2)) 
(15) is not antisymmetric. Eq. 13 describes the superpo¬ 
sition state of system electron being in upper and lower 
potential wells, while entangled with its spin degree of 
freedom. When 9m = 0 or tt, there is no splitting of 
the wave packet and eq. 13 reduces to the correct form 
describing the unsplit wave packet inside the potential 
well. 


Inside the double well potential (excluding the bound¬ 
ary points, as derivative of the wave function does not 

2 2 

exist there) Hqs = and it is easy to 

check, using ^0[z — zq) = S(z — zq), that the state 
|4>os„(z)) given in eq.l3 is an eigenket of Hqs with eigen- 

value En = ■ 

We can make the system electron to enter the eigen- 
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state |$osn(2)) of Hqs, by choosing = ^ = 

^ilBi^T)"^ and following an argument similar to that 


carried out from eq. 28 to 31 with the limit Vq —t oo (here 
we also take the limit Ap^ —t 0. However, this is not nec¬ 
essary if we follow a different technique as explained in 
para next to next to eq. 31). Hence, in principle we can 
push the system electron into one of the eigenstates of 
Hqs- 

Wavefunction in eq. 13 is Square Integrable unlike the 
plane waves exp(^pzz) and hence represents the actual 
physical state of system electron. Hence there is no 
need to construct a wave packet using the eigenkets 
{|‘l’osn(z))}, to represent the physical state of system 
electron inside the potential well, unlike the situation 
when the electron was free. This justifies our assump¬ 
tion that system electron enters the eigenstate given in 
eq. 13, inside the potential well. This can be further 
justified by the fact that, an electron in a bound state of 
hydrogen atom can exist in one of the eigenstates, even 
though it was in a superposition of infinite number of 
plane waves exp(|p.f^ when it was unbounded. 

Completeness of eigenstates of Hqs: Using eq.l3, 
taking the inner product of |4>osfc(.z)) with |4’osn(z)) 
which involves integration w.r.t z from —oo to -boo, we 
get ($osfc(z)|4)os„(z)) = Skn, where 4n is the kronecker 
delta function, and hence orthonormal. In eq. 13 if 
we introduce a relative phase (e*’^) between two terms 
(which is equivalent to changing 9m to —6m), we get a 
normalized solution which is linearly independent from 
|4>osn('2))- Then using Grahm-Schmidt orthogonaliza- 
tion, the solution becomes: 
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|$Oan{^:)) =sin^| 0 ) e{z - {b - - 0{z - {b + -)) \ - sin {koanib + - - z)) 


-cos^e*'^’"|l) 9{z + {b +-)) - 6{z + {b --)) \ - sin {koan{b + - + z)) 


(15) 


where kgan = n = 1,2,3,.... One can verify 
that |$oon(-z)) is an eigenket of Hqs, with eigenvalue 
same as that corresponding to the state |$os„( 2 ;)), and 
hence degenerate with it. One can also verify that 


(^0afc(2:)|^’0an(2:)) = 4n and (4*08* (^) | $0an ( 2 :)) = 0 for 
all k,n and hence mutually orthogonal. Using the rela¬ 
tion: 


^l^sin ( —z) ^^sin ( —z') [6{z) — 6{z — a)] [0{z') — 9{z' — a)] = 5[z' — z) [9{z) — 9{z — a)] [9{z') — 9{z' — a)] 

n—\ 

(16) 


where, (5(z' —z) is the Dirac delta function, one can verify that: 


E ( l^ 08 „(^))(<i> 08 n( 4 )| -f \^ 0 aniz)){^ 0 an{z')\j = |0)(0|4^' - ^) [o^Z - (b - ^)) - 9{z - (b + ^)) 
Oiz' - (6 - ^)) - 9iz^ - (6 + ^))] + |l)(l|4z' - z) [9{z + (b + ^)) - 9{z + {b - ^)) 


9{z' + {b+^))-9{z' + {b-^)) 


dz' Y.{\^ Osn (Z))(4>08n(^')l + l‘&0an(^))($0an(^')l ) = 1 (17) 


where, 1 is 2 x 2 identity matrix, there by satisfying the 
completeness relation. 

To measure a given system protectively, system should 
be in a non-degenerate eigenstate of system Hamilto¬ 
nian. But the state |4>os„(z)), which we want to mea¬ 
sure protectively, is a degenerate eigenstate of i?os, as 
discussed in previous paragraph. If we perturb the sys¬ 
tem by applying a small negative potential in the region 
— {b—^)<z<{b—^), degeneracy will be lifted due to 
tunneling. Adding a small negative potential to a very 
large but finite potential(which is infinite for all practical 
purposes), lowers the net potential slightly, there by al¬ 
lowing tunneling. Instead of solving perturbatively (for 
perturbative treatment see [ 10 ]), we are going to treat 
the potential in the region — (6 — |) < z < (6 — |) to 
be Vo(> 0 ), which is finite and constant, and solve ex¬ 
actly. If necessary, later we can take Vq to be very large 
but finite, to satisfy the condition of small perturbation 
to Vi{z) (11). Following treatment is similar to that on 


pp460-62 of [10]. We are going to find the eigenstates 
\4‘sn{z)), \(j)an{z)) of new System Hamiltonian: 

Hs = §^+VAz) ( 18 ) 

with eigenvalues Egn and Ean respectively, which are less 
than Vq. Vf{z) is same as V7(z) given in eq. 11 but 
Vf{z) = Vo in the region —{b— §) < z < (6 — |) instead 
of oo. 

As the system electron is in superposition of being in 
both wells, it can tunnel from both wells into the region 
— {b — a/2) < z < {b — a/2.). As the split wave packets 
tunnel into the region —(5 — a/2) < z < {b — a/2), they 
get united and hence disentangles from spin d.o.f. This 
is justified by the fact that, if we reverse the direction 
of inhomogeneous magnetic field (see above eq. 49), the 
split wave packets start moving towards the origin, get 
united and hence disentangles from spin d.o.f. With this 
requirement, |4>osn(2:)) (13) suggests the following form 
for one of the eigenstates of Els with energy Eg < Vq: 
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|(;i«(z)) = cos^lO) 6 »(z-( 6 -“))- 6 l(z-( 6 +“)) A sin{ks{b +- z)) 
+ sin^e*'^"‘|l) 6 l(z+( 6 +“))- 6 »(z +(5- “)) A' sin {ks{b + ^ + z)) 
+ (cos %|0) + sin |1)) \e(z + (6 - “)) - O^z - (b - “))! [B + B' 


(19) 


where, qs = ^J^iVo — Eg) = — k^ > 0. We have 

related energy Eg of the state with its wave number kg 
by the relation: 


Eg 


2M 


( 20 ) 


Exponentially decaying wave corresponds to tun¬ 

nelling from upper potential well (i.e., potential well on 
positive z-axis) into the region —{b— |) < z < (&— |) and 
g-qs{z-{-b)) corresponds to that from lower well. Fur¬ 
ther, |<i)osn(z)) (13) suggests to take A = A' = Ag and 
B = i?' = in eq. 19. Evanescent wave reduces to 
the form: B -|- B' = Bg cosh(gsz), 

where Bg = 2B'ge~^‘‘^. It is evident that in the limit Vq 
(hence Qg) going to infinity, both evanescent waves vanish 


and we recover the state in eq. 13 as required. Demand¬ 
ing that the spatial wave function and its first derivative 
must be continuous at z = 6 — |, we obtain the following 
constraint equation: 

tan(fcsa) =- / ^ coth(v/a^ - k'^{b- “)) ( 21 ) 

-kl 2 

This follows from the fact that Ag and Bg cannot vanish 
simultaneously, else we get trivial solution |(^s(z)) = 0 
everywhere. Also in obtaining eq. 21 we dropped unit 
step function which is justifiable as it is used just for 
the sake of convenience. Similar joining conditions at 
z = —{b — |), also gives same constraint eq. 21 . 

Now, |d)oan(- 2 )) (15) suggests the following form for an¬ 
other eigenstate of Hs (which is linearly independent 
from \(j)g{z)) 19) with energy Ea < Vq- 


|(()a( 2 :)) = sin^lO) 6 »(z-( 6 - “)) - 6 l(z-( 6 -b “)) A sin {ka{b + ^ - z)) 
-bcos^e*'^"‘|l) e{z + {b+^))-e{z + {b-^)) A' sin {kaib + ^ + z)) 
+ (sin %|0) + cos |1)) [0(z + (6 - “)) - 0(z - (6 - “))! (B + B' 


( 22 ) 


where, Qa = Ea) = ^- k^ > 0, = 

h'^k'^/(2M). |$oan(-z)) (15) suggests to take A' = —A = 
—Aa and B' = —B = —B' . With this, evanescent 
wave takes the form: B + B' = 

Ba sinh(ga-z )5 where, Ba = 2 il(,e“‘^“*’. In the limit 
Vq —t 00 , \(j)a{z)) (22) reduces to |d)oan(z)) (15) as re¬ 
quired. Applying joining conditions at z = 6 — | in a 
manner similar to previous case, we obtain the following 
constraint equation: 

tan(A:aa) =- , tanh - kl{b - ^)) (23) 

Joining conditions at z = —{b— |) also leads to same 
constraint eq. 23. Equtions 21 and 23 can be solved 
graphically for kg and ka respectively, to obtain kgn and 
kan as roots, there by quantizing the energy levels. 


As constraint eq.s 21 and 23 are not identical in form, 
kgn will not be equal to kan- Then, using relation 20 and 
Ea = h‘^ka/{2,M), we obtain Egn ^ Eam there by lifting 
degeneracy. In the limit Vq —t oo, we obtain from eq.s 21 
and 23: kg^a ^ as required. As Qgn.an > 0 we have 
Egn,an < Vq- Eigenstates of Hs, given by eqs. 19 and 
22 , take the following form after quantization: 


\ 4 >sn{z)) = \4>0sn{z)) J- \m)Bgn COSh(gs„z) 

o{z + ib-l))-e{z-{b-^)) 

If^aniz)) = \(l)0aniz)) + (sin^|0) -b COS ^11)) 

Ban sinh{qanz) &{z + {b - ^)) - diz - {b - ‘^)) (24) 


where, |m) is given by eq. 8 , and: 
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\^Osn{z)) = COS ^|0) e{z - {b- “)) - 9{z - (6 + “)) A^n sin + “ - z)) 


+ sin 11) [eiz + (6 + “)) - 0(z + (6 - “)) 


Asn sin {ksn{b +^+z)), 


l^ 0 an( 2 :)) = sin^lO) 6»(z- (5- “)) - 6»(z - (&+ “)) Aa„ sin {kan{b + ^ - z)) 


id 

— COS —e ^ 

2 


^11) + (6+“))-0(z+ (&-“)) 


^an Sin(fcari(&- 


^)) 


(25) 


In the limit Vq —)■ oo, we obtain \(psniz)) —)■ |$osn( 2 :)) 
and \4>an{z)) —t |$oan(-z)) as required. One can ver¬ 
ify that {4>ak(.z)\4>sn{z)) = 0, for all k,n. However, 
(^sfc(-z)|(/'sn(z)) fy 0 and ((('ofc(-z)|(/>an(-2:)) fy 0, for fc fy n. 
We can make them zero via Grahm-Schmidt orthogo- 
nalization. Arbitrary constants Agn, Aan, Bgn and Ban 
can be fixed by the requirement that each of the states 
{\(j)sn{z)),\4>ak{z))}, for all s, fc, be normalised and sat¬ 
isfy completeness relation along with the states {\4>'i{z))}, 

where \4>iiz)) = \i^)cos{i^^z) is an eigenstate of 

Hs with eigenvalue E[ > Vq. Also {(j)[{z)\(j)sn,an{z)) fy 
0. Taking the state of our main interest, fysi(z)), 
as reference, we can orthogonalize the complete set 
{\(j)sk{z)), \(j)ak{z)),\(l)[{z))} via Grahm-Schmidt orthogo- 
nalization. The set {\(j)sk{z)),\(l)ak{z)),\(l>'i{z))}, after or¬ 
thonormalization, should satisfy the completeness rela¬ 
tion analogous to that in eq. 17, i.e.. 


r M 


J dz' ^ \^\(j)gkiz)){ 4 'sk{z')\ -I- \(l)ak{z)){(l)ak{z')\j + 

-oo ^=1 

OO 

oo 1 „ oo 

^ mz))mz')\ = / dz'^fy"(^))(</''.V)l = l(26) 


/=2M+1 


where, fysfc(z)) = |(/> 2 fe-i(z)), |<(>afe(fy) = |(/' 2 fe{^)), 

\^[{z)) = |</.r(z)) (27) 

We can justify this by the fact that, in the limit Vq —)■ oo 
we recover the completeness relation in eq. 17, which cor¬ 
responds to Vi{z) (11). By perturbing Vj with a small 
negative potential, we obtain the case of finite potential, 
Vf{z) (18). As perturbation doesnot destroy the com¬ 
pleteness property of the set of states being perturbed, 
completeness relation in eq. 26 is justified. 

Protective Measurement: An important require¬ 
ment to do protective measurement is, system should be 
initially in a nondegenerate eigenstate of system Hamilto¬ 
nian, Hs- Energy of system electron, just before trapping 
it in double well potential is: 

where, poz is the average momentum, Ap^ is the standard 


deviation in momentum and 

AA = ^{^pI ± ‘^iBi^rAp,) (29) 

is the uncertainty in kinetic energy of system electron. 
Let us consider the ground state \(j)si{z)) (24) (as Egn 

is less than Ean [10]). From eq. 20, Egi = But 

eq. 21 says, depends on the potential well parame¬ 
ters a, 6, Vb- To push the system electron into the state 
|(()si(z)), following condition must be satisfied, 

Egi{a,b,Vo) = (30) 

Even if the values of a, 6, Vq in Esi{a, b, Vq) are fixed by 
some requirements like: compensating increase in Az as 
Apz is squeezed (see below eq. 31), etc., constraint eq. 
30 can still be satisfied by varying the free parameters Bi 
and r. Hence, we can say, energy of system electron, just 
before trapping is Egi + AE. In general, state of system 
electron after trapping it in double well potential will be, 

|(/>d»p) = Ci|(/)^i(z)) +^Cify7(z)) (31) 

where, \(j)”{z)) is given by eq. 27. However, system elec¬ 
tron starts with zero average momentum as mentioned 
below eq. 9. At t = 0, AzAp^ = fi./2. In principle 
we can start at t = 0 with squeezed wave packet such 
that Apz —>■ 0. Switch on the inhomogeneous magnetic 
field for an interval of time r such that it gains an av¬ 
erage kinetic energy Esi{a,b — a/2 -|- e, Vq) (e > 0, to 
avoid blowing up of cot hyperbolic function in eq. 21) 
as given by the constraint eq. 30. Now adiabatically 
change the Hamiltonian from p^/(2M) to Hs (18) with 
b = a/2 + e. By this, energy of system electron will 
not change, as we are introducing potential barriers only. 
Width ‘a’ of well is suitably chosen so as to take care of 
increase in Az due to squeezing of Ap^.. In principle ‘a’ 
can be arbitrarily large. With this we obtain: AE 0, 

I Cl I —>■ 1 and jcij —>■ 0 (z fy 1) with state of system elec¬ 
tron being |^si(z)) corresponding to 6 = a/2 + e [14]. 
Now again adiabatically change the Hamiltonian from 
Hs with b = a/2 -|- e to Hs with b = a/2 -|- e -I- Al, where 
Ai{> 0) is large but finite i.e., separate the two wells 
very slowly. Hence, in principle, ‘6’ can also be made ar¬ 
bitrarily large. Here, energy of system electron changes. 
Hence, by adiabatic principle (see below eq. 36), system 
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electron will end up in the state \(l)si{z)) corresponding to 
b = o/2+e+A^. Hence, in principle, we can push the sys¬ 
tem electron into the nondegenerate eigenstate, |0si(z)), 
of system Hamiltonian Hs- 

Or, it is sufficient for the system electron to have av¬ 
erage kinetic energy (which is 2 ® (7'®*!''’)^) sufficiently 
greater than the ground state energy Egi, which is ex¬ 
actly known. As, all systems have natural tendency to 
settle down to their ground state, as it is most stable and 
has least possible energy, if we wait sufficiently long, sys¬ 
tem electron will settle down to the ground state \(j)si{z)) 
[15] [16]. Hence, in this method we need not take the 
limit Apz —>■ 0, unlike in the previous method. Hence, 
by either of these two methods, we can push the system 
electron into a nondegenerate eigenstate |0si(z)) (24) of 
system Hamiltonian Hs (18). Even if we trap system 
electron in the eigenstate \4>sn{z)), n > 1, by satisfy¬ 
ing the constriant: E^n = , it won’t be sta¬ 

ble due to vaccum fluctuations and ultimately decays to 
ground state, \4>si{z:))- 

We now carry out protective measurement on system 
electron trapped in the double well potential, using a 
probe electron which passes in-between the two wells. 
As it passes, it interacts with system electron through 
coulomb potential energy. Even though we have trapped 
system electron in 1-D potential well, it can still inter¬ 
act with probe electron moving in x-z plane(let positive 
x-axis be along positive V/-axis in Fig. 1), via coulomb 
potential energy. Trapped electron is analogous to line 
charge density (Coulomb per meter), where the potential 
well is opaque to matter wave(hence charges cannot es¬ 
cape from metal surface) but transparent to electric field 
and hence they can interact with other free charges out¬ 
side the well. Only trapped system electron sees infinite 
potential barrier in between the two wells, but not the 
probe electron [17]. 

As stated in [9], adiabatic condition required to do pro¬ 
tective measurement will be satisfied, assuming system 
electron is initially in the nondegenerate energy eigen¬ 
state \(j)sn{z)), provided the probe electron crosses the 
potential well in a time T, which is large compared to 
h/ \Esn — Ean\ where, Egn, Ean are the energies of the 
states \(j)sn{z)}, lifaniz)) respectively. As the probe elec¬ 
tron is moving very slowly, we can neglect the magnetic 
field induced due to motion of charge. In a more rigor¬ 
ous calculation we can treat it as low energy Quantum 
Electrodynamical (QED) problem. However, here we are 
going to approximate it as quantum electrostatic prob¬ 
lem. Let the probe electron have an initial constant mo¬ 
mentum, PaOxi, where i is an unit vector along x-axis. 
As it measures the system electron protectively, it gains 


component of momentum along z-axis, and hence drifts 
along, say, positive z-axis. As a result, potential energy 
of interaction, H°■„^ = changes, where f 

is the position vector of system electron and ra{t) that 
of probe electron. However, if the separation between 
wells, 2(6 — a/2), is sufficiently large (as justified be¬ 
low eq. 31) and time interval T relatively small (i.e., 
T >> h/ \Esi — Eai\ still holds) [18], we can neglect the 
change in Hf^^ during the time interval T and approxi¬ 
mate it as: 


HI 


1 

47reo \z - Za\ 


(32) 


with, center of wave packet of probe electron: Zao = 0. 
It is as if probe electron is fixed exactly in between the 
two potential wells. 

As the separation between two wells, 2(6 — a/2), can 
be made arbitrarily large as justified below eq. 31, 
1/ |z — Za[ can be made arbitrarily small for every z such 
that: (6 — a/2) < z < (6-1- a/2) and —(6 -|- a/2) < z < 
— (6 — a/2). Hence H/,^^ can be made arbitrarily small 
compared to \Esn — Ean\ (also see last para of [19], and 
second para of [20]), there by satisfying the weak interac¬ 
tion requirement for protective measurement, i.e., is 
so weak that it cannot cause transition between the states 
{(/sn) and l^an). We can vary \Esn - Ean\ by varying Vq. 


As we are interested only in the change in z-component 
of probe electron’s momentum, we can neglect its other 
components. With this and the electroststic approxima¬ 
tion in eq. 32, we can take free evolution Hamiltonian of 
probe electron, Ha =p^^/(2M) = 0 [21]. Hence, System 
and probe electrons evolve under the Hamiltonian: 


H = Hs + 


1 

47reo [z - Za[ 


(33) 


where, Hs is the system Hamiltonian given in eq. 18, z is 
the position coordinate corresponding to system electron 
and Za that of probe electron. State of combined system 
at time t = T is: 

|«T)) = exp (^(Hs + 

(34) 


where, \(j)si{z)) is given by eq. 24, Q{paz — 0) is the 
Gaussian wave packet of probe electron in momentum 
space centered at PaOz = 0. Eollowing Hari Dass and 
Tabish [15], we insert completeness relation similar to 
that in eq. 26, to obtain: 


lan 



00 



00 

dz' exp 



1 

47reo |z - Za\ 


T ]5{Za){\^k{z)){^k{.z')\)\4>sl{z'))Q{paz) 


(35) 


— 00 


— 00 








where, hat on z, Za stands for operator, S{za) is the Dirac 
delta function and \(j)si{z)) has been changed to \(j)si{z')) 
because see[22]. If, A\ai) = ai\ai) then f{A)\ai) = 
f{ai)\ai). Using this property we can push S{za) to the 
extreme left. Then we obtain: 

(hs + - 1 ) l^fc(^)) = Ek{zaMk{z)), where, 

\ iTT€o\z-Za\J 

Ek{za) = ^ \Mz)) (36) 

As discussed in the beginning of ‘Protective Measure¬ 
ment’, criteria of adiabaticity and weak interaction be¬ 
tween system and probe electrons are satisfied. Accord¬ 
ing to Adiabatic theorem, at t = 0 if the system electron 
is in a nondegenerate eigenstate of the system Hamilto¬ 
nian Hs (which we take to be the ground state, |(/)si(z))), 
then, if Hs changes very slowly to H = Hs + at 
t = T, such that T » (hj — A'ail) (adiabatic ap¬ 
proximation, eq. 10.15 of [23]), then \(j)si{z)) changes to 
the corresponding nondegenerate eigenstate of H (with 
Za replaced by Za, as in eq. 36) i.e., ground state |0i(0)) 
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(pp327-29 of [23]). This is true for arbitrary Hf^^ (see the 
generalization on p330 of [23] ). Further, as explained be¬ 
low eq. 32, if the separation between wells is arbitrarily 
large, then, Hf^^ will be arbitrarily small (also see [20]). 
Then, by nondegenerate perturbation theory, correction 
to all orders to the state |0si(z)), due to perturbation 
Hf^^, goes to zero [24] and we obtain: |<^i( 2 :)) tends to 
\(j)si(.z)). Generalizing to an arbitrary energy eigenstate, 
we obtain: \(j)k{z)) tends to |0).'(z)), where \4>k{z)) is de¬ 
fined in eq. 27. Using orthonormality of states, i.e.. 


dz'{c^'l{z')\cl^Az')) = 

OO 

/ dz'{(l)'l{z')\(f)ak{z')) = 
dz'{<^'l{z'M{z')) = 


1, if i = 2fc - l,fc < M 
0, otherwise 


— OO 
OO 


1, if i = 2k, k < M 
0, otherwise 

1, Ai = l,l>2M + 1 
0, otherwise 


(37) 


and using eq.36, |C(2^)) in CQ- 35 reduces to. 


\({T)) = j dza5{za)eyiY>{-^{Hs)4,^^T)eyi^(^-^^^^ 

— OO 

_ 2 , 

= exp ( —(i7s)0^jr) |(;isi(z))exp 


Let, /(za) = ((/>si(^)|]— ^^Msiiz)) (39) 

\Z - Za\ 

Using the definition of wavepacket (pl462 of [10]), we 
obtain: 

exp ^ 

OO 

^ + (40) 

— OO 

Wave packet of probe electron in position space is cen¬ 
tered at Zao = 0 (see below eq. 32) and center won’t 
change with time, as PaOz = 0 (see above eq. 10). Ex¬ 
panding f(za) around Za = Zao = 0, we obtain, f{za) = 
f{0) + f'{0)Za + Oizl), where, /'(O) = -£^f{za) U„=o- If 
the spreading of wave packet is small during the time in¬ 
terval T [25] [9], Q{za), a Gaussian wave packet, is appre¬ 
ciable only in the interval: —e < Za < e, e > 0, such that 
and higher order terms are negligible i.e., 0(e + 5) ~ 0, 
where, 5 > 0. Hence in integral (40), f{za) can be ap¬ 
proximated as: f{za) ~ /(O) -I- f'{0)za- Eq. 40 reduces 


(</>sl(^)l)- ^ - Msl{z))T\\cj)si{z))e{paz) 

IZ-Zal J 

(<(>.1 (^) I 0 ( Pa .) ( 38 ) 

ft 47reo Iz - Zal ) 


to [26]: 

Finally we obtain, 

|C(r)) = exp {El (^Hs)^^, + ^eV(0))r) 

|<(>.i(^))0(^Pa. + ^eV'(O)r) (42) 

State of system electron is unaltered while probe electron 
has gained an average momentum: ~ 4 ;^e^/'( 0 )T, and 
we are now going to show that f(A) is proportional to 
cos0m- Substituting eq.24 into eq.39, we obtain: 

f{Za) = (</>0sl(2:)|)- - - -\(j)osl{z)) 

\z- Za\ 

+ [ dz = h{za) + g{za) (43) 

J ^a\ 
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In eq.s 43 and 45, integration variable is z, and hence Za 
is a constant as far as integration is concerned. We want 
to know the rate at which the function f{za) changes, 
w.r.t Za, at Za = 0. Hence w.r.t Za = 0 (the value of our 
interest) we can write: 


\Z - Za\ 


Z — Za, for z > 0 
— (z — Za), for z < 0 


(44) 


However, we want /'(O) but not /(O). Hence, we are 
going to first differentiate f{za) w.r.t Za and then put 
Za = 0. Using (44) and substituting for \(j)osi{z)) from 
eq.25 and integrating w.r.t z from —oo to +oo, we obtain, 


h{Za) = ((()0sl(2)| I -. |(;iosl(^)) 

\Z - Za 


b-\- 7 


= Ai COS -Y J 

b+f 

sin^ ^ 

b-^ 


sm‘^{ksi{b+ § - z)) 


sin^(fcsi(&+ ^ - z)) 


Z + Za 


(45) 


Differentiating h{za) w.r.t Za in eq. 45, using Leibnitz 
rule for differentiation under the integral sign, and then 
taking Za = 0, we obtain. 


b+f 


/i'(0) = cos 9m J dz 

b-% 


sin^(fcsi(6+ I - z)) 


(46) 


Similarly, 


b-§ 


9i.Za) = 


/ cosh^((7siz) 


z + Za 


cosh^(g^iz) V^^^ 
- -Za / 


To avoid pole at z = Za = 0, we have introduced the 
limit. As before, differentiating g{za) in eq. 47, w.r.t Zq, 
and then taking Za = 0, we get ^'(0) = 0. .’. using eq. 43, 
we obtain /'(O) = h!{Q) (46). From eq. 42 we obtain the 
final average momentum of probe electron (whose initial 
average momentum was zero) to be: 


final 

PaOz = 


47ren 


b-\- % 


COS0„ 




b-% 


In principle we can measure momentum of probe electron 
(and hence dm) with arbitrary precision [19]. We note 
that, both T —>■ oo and separation between two poten¬ 
tial wells tending to infinity are necessary. In the latter 
limit, the integral in eq. 48 tends to zero (as integration 
is over only upper well), there by keeping the product, 
T X integral, finite. We also note that, in an exact cal¬ 
culation (i.e., without doing any approximations), the 
functional form of eq. 48 may change, but still it will be 


a function of 9m- In eq. 48, except cos 9m all other terms 
are greater than zero. If the unknown 9m happens to be 
zero, there is no splitting of wave packet and the system 
electron will be in upper well for Bi < 0, according to eq. 
10. Hence, probe electron experiences maximum down¬ 
ward force(repulsion) and hence gains maximum negative 
momentum, which is consistent with in eq. 48. If 

= f, wave packet splits equally and hence probe elec¬ 
tron experiences equal and opposite forces, gaining zero 
net momentum. in eq. 48, also gives same result. 

Finally if 9m = tz, again no splitting of wave packet, but 
this time system electron will be in lower well. As a re¬ 
sult, probe electron experiences maximum upward repul¬ 
sive force, gaining maximum positive momentum. Again 
in eq. 48 is consistent with it. 


If we just want to discriminate between two nonorthog- 
onal states, say |0) and (|0) -|- |I))/-\/2, then we have suc¬ 
ceeded in discriminating, as it requires only the knowl¬ 
edge of polar angle 9m- However, if we want to clone 
the unknown state completely, we need to still findout 
its azimuthal angle, (pm (8). As the state of system 
electron is unaltered at the end of protective measure¬ 
ment (42), and because the trapping process was adi¬ 
abatic (see below eq. 31), it must be reversible, and 
hence we can get back the state in eq. 10, by adiabat- 
ically untrapping (procedure for untrapping is same as 
for trapping (below eq. 31), but in reverse). Now we 
can recombine the split wave packets [20] by reversing 
the direction of inhomogeneous magnetic field which was 
used initially for splitting i.e., subjecting it to a field: 
B = BB = —BiX i {—Bq Biz)k (refer [11]), where, 
Bi > 0 (note that while splitting, Bi was < 0). We re¬ 
cover the initial state, \tp{0,Pz)) = \'in)x{Pz — 0) (10), 
which had Ap^ —)■ 0, but with AzAp^ > h/2 [27]. Polar¬ 
ization, m, of unknown spin state is unaltered w.r.t lab 
frame. Previously we had applied inhomogeneous mag¬ 
netic field along positive z-axis of lab frame. Now we 
apply it along a unit vector, 

h = sin 9n cos (pni + sin sin <pnj cos (49) 

where, 9n ^ Q (zero is the value corresponding to previ¬ 
ously applied direction of inhomogeneous magnetic field) 
and (pn 7 ^ 0. From eq.s 49 and 9 we obtain the following 
constraint equation: 

rh-fl = cos 9mn = cos 9m cos 9n-\-COs{(pm — (pn) siu 9m sin 

(50) 

Component of unknown spin angular momentum opera¬ 
tor along n is, Sn = S.n, where S = S^i + Syj S^k = 
Smizi, Sm is the component along itself. Eikenkets of 
operator Sn are: 

l-k)n = cosyjO) -f sinye*'^"]!) 
and ]-)„ =-sinylO)-l-cosye*'^"]!) (51) 
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with eigenvalues | and respectively. In the eigenbasis 
of Sn, unknown spin state (8) can be decomposed as: 

|m) =cos^|+)„ + sin^e*'^’"'‘|-)„ (52) 

where, 0mn is given by eq. 50 and (pmn is the az¬ 
imuthal angle of the unit vector m, in a coordinate sys¬ 
tem rotated w.r.t the fixed lab frame, in which n is 
along it’s positive z-axis. Inhomogeneous magnetic field 
along h is switched on for an interval of time r. Sys¬ 
tem electron evolves under the interaction Hamiltonian 
Hint = -fl-Biqnfi = -jBiqnSn, where is the dis¬ 
tance measured from origin along the direction n. State 
of system electron after interacting with inhomogeneous 
magnetic field is given by: 

|'0n('r,p„)) = cos ^|-f)„x(p„ - ibJ^t) 

-ksin^e*'^’""|-)„x(p„-f (53) 

where, Pn is the component along n of momentum of sys¬ 
tem electron i.e., momentum conjugate to qn- As before, 
trapping the system electron in a double well potential 
and measuring it protectively with a probe electron, we 
obtain the value of cos 0mn ■ Hence we obtain the con¬ 
straint eq. 50 with the only unknown (j)Tn- However we 
cannot obtain the value of (j>rn unambiguously from a sin¬ 
gle constraint equation for the following reason: Consider 
the constraint eqation, cos?] = c, —l<c<l. If?] takes 
a value in the interval [0,7r] (just like 0m)^ then there 
is only one value of ?] which satisfies cos ?] = c. On the 
other hand, if ?] takes a value in the interval [0, 27r) (just 
like then there are two values of rj which satisfies 
cos rj = c. Hence we require one more constraint equa¬ 
tion independent of the constraint eq. 50, to find out 
(l)m unambiguously. This is justified by the fact that, to 
unambiguously characterize an arbitrary orientation in 
space, we require three Euler angles, which are linearly 


independent of each other. Whereas, here we are trying 
to characterize the unknown direction, m, with just two 
parameters 0m and (ftm, and hence the ambiguity. 

After obtaining the value of cos0mn, we shall again 
recombine the split wave packets, as before. Applying an 
inhomogeneous magnetic field along the direction I given 
by, 

I = sin 01 cos (l>ii + sin 0i sin -|- cos Oik (54) 

where, 0i ^ 0,(0, 0n are values corresponding to previ¬ 
ously applied directions of inhomogeneous magnetic field) 
and (j}i 7 ^ 0, By this we obtain another independent 
constraint equation: 

m.l = cos 0ml = cos 0m COS 01 + COs{(j)m — </'?) sin 0m sin 01 

(55) 

Through protective measurement we can obtain the value 
of cos 0ml- Solving eq.s 50 and 55 for tpm and picking 
out the solution common to both of them, we obtain an 
unambiguous value of (pm- 


III. MISCELLANEOUS 

Discriminating between |0) and |-l-)(= (|0) -I- 
|l))/v^) even when the initial state is in a lin¬ 
ear superposition of nondegenerate eigenstates of 
Hs (18): Let the initial state of system and probe elec- 

tron be: |C'(0)) = £ Cj|(('"(z))0(paz), where, \(p'l{z)) is 

i=l 

defined in eq. 27. Hence we need not worry about push¬ 
ing the system electron into one of the eigenstates of iLg, 
unlike in cloning. Hence discrimination is more feasible 
than cloning. Allowing |C^(0)) to evolve under the Hamil¬ 
tonian H given in eq. 33, for an interval of time T, and 
inserting completeness relations we obtain: 


\C{T)) = ^^0^ J dza J dz'^exp (^-j^i^Hs + ^^ . ^ ^T^d{za)\pkiz)){pkiz')\\(p'Iiz'))‘d{Paz) (56) 

^ — ^ -OO -OO ^ = 1 


This is similar to the state in (35) but with an extra 
summation over the index As we are doing Protec¬ 
tive Measurement, criteria of adiabaticity and weak in¬ 
teraction are satisfied and hence as before (below eq. 36) 
we can take, \(pk{z)) ~ \(p'p.{z)). Proceeding in a fashion 


similar to that from eq. 35 to 48, we obtain, 

IC'(r)) = £ c, exp T) 

|</>"(z))0(^Pa. + ^eV'(O)r^ (57) 
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where, 


f)+§ 


fi{0) = cos9m Aj J dz 


sin^{ki{b + I - z)) 


(58) 


b- § 


If the given state is |+), then 9m = ^ and hence /((O) = 0 
for all As a result, there is no entanglement between 
system electron and probe electron (eq. 57) and probe 
electron gains zero momentum. This makes sense, be¬ 
cause, when 9m = ^, wave packet of system electron 
splits equally. As a result probe electron experiences 
equal and opposite repulsive forces, gaining no net mo¬ 
mentum. However if the state is |0), then 9m = 0 and 
hence /((0) is strictly greater than zero as being nor¬ 
malization and completeness satisfying constant, cannot 
be zero and ki (20) which corresponds to the energy of the 
level, cannot also be zero. Hence in this case, probe 
electron gains non-zero momentum (of course, in this case 
there will be collapse upon measurement as there is en¬ 
tanglement, but it does not matter as we are interested 
only in discrimination). Also, wave packet of system elec¬ 
tron doesnot split and hence it will be in the upper well. 
As a result probe electron experiences downward repul¬ 
sive force, there by gaining downward momentum. By 
measuring the momentum of probe electron (which can 
be measured as precisely as we want, at the cost of loos¬ 
ing information about its position, which is not required 
for discrimination, any way) we can always discriminate 
between the states |0) and |-|-) unambiguously. However, 
we obtained eq. 58 under the approximation given in 
expression 32 and a few other approximations. Suppose 
there is a small correction A/'(0) to /((O) in eq. 58, un¬ 
der an exact calculation. By choosing a and b suitably, 
we can make the integral along with Af in eq. 58 as 
large as we want. When 9m = f, /((O) = A/((0), where 
as when 9m = 0, /((O) >> A/((0). Hence, we can still 
discriminate between |0) and |-|-) always. However, we 
note that, A/'(0) may also increase as the integral in 
eq. 58 increases, in which case above argument breaks 
down. Also it may happen that under an exact calcula¬ 
tion, there might be entanglement between system elec¬ 
tron and probe electron even for 9m = f > in which case 
we cannot discriminate, as probe electron also gains non¬ 
zero momentum. One has to do exact calculation and 
see. 


How Not to Clone: It’s not possible to clone an ar¬ 
bitrary unknown spin state |rTi) (8) by protecting it with 
a strong homogeneous magnetic field (Bq). As the di¬ 
rection, m, of unknown spin polarization is not known a 
priori, we cannot apply Bq along m, in general. Hence we 
choose the direction of Bq to be along k (positive z-axis 
of fixed lab frame) i.e., Bq = Hq/c. Apply an inhomo¬ 
geneous magnetic field ^qnfi, where, is the distance 
from origin measured along h (49) , ^ is the field gradient 
in Tesla m~^ and in this context, T is a dimensionless 
number which controls the gradient strength. As, Bq can 


be arbitrarily large, we can neglect pl/{2M) compared to 
I — /i.i3gfc|. Hence, |m) evolves under the Hamiltonian, 

H = -fl.BQk - fl-B—n = -fl-iBok + B~n) 

= —"fS.B = —jBSb (59) 

where, S is the total spin angular momentum of unknown 
spin. Sb = S.B, where H is a unit vector along resultant 
magnetic field B. In the eigenbasis of the operator Sb, 
state of unknown spin-^ particle at time t = T is, 

\M'{T)) = cos^e"5'^®|-b)Bexp(-^(-7H)^T)x(pn) 

-t sin |-)b exp(^(- 7 H) ^T)x{Pu) (60) 

where Pn is the componenet of momentum of unknown 
spin-i particle along h, cos 9Bm = B.m. In the limit 
T —)■ 00 , which corresponds to adiabatic and weak in¬ 
teraction limit, we obtain: B = \B\ ^ Bq + ^(7„cos0„. 
But cos9n = |(S'„)|o), where, Sn = S.fi and |0) is the 
eigenket of Sz- In this limit, state \M'{T)) becomes: 

|M'(T^ oo)) = 

g-5(0i3-7So7’) cos ^^\+)Bx(^Pn “ 7H,(S'„)|o)^ 

sio^^\-)gx(pn + lB,{Sn)\0)^ (61) 

In the limit T —>■ oo, 


(62) 


where, cos 9mn is given by eq. 50. Consider the case 
where, the direction of unknown spin polarisation, m, 
happens to accidentally coincide with that of static field 
direction, k, and hence ‘rightly protected’ (this is the case 
considered in [9] as mentioned by themselves in [28]). In 
this case 9m = 0 and we obtain cos 9Bm ~ I ^ 9Bm ^ 0, 

B ~ m = k, |-|-)b ~ |m) = |0) and the state 61 becomes, 

\M'{T ^ oo)) = - 7i3*(5n)|o)) (63) 


There is no entanglement in eq.63 and the unknown 
spin state is unaltered. Upon measuring the momentum 
of spin-i particle, we obtain: 7 i?i(S'„)|o) = xBi{Sn)\m)- 
Consider now the case where, the direction of unknown 
spin polarisation, m, donot comcide(which is highly 
highly probable) with that of static field direction, k, and 
hence ‘wrongly protected’. In this case, cos9Bm ^ I as 
evident from (62) (this is true even if 9mn = 0 i.e., to ac¬ 
cidentally coincides with n ) and hence spin and spatial 
d.o.f are entangled as given by eq. 61. Upon measur¬ 
ing the momentum of particle in momentum basis, we 
obtain either, momentum xBi{Sn)\o) and spin state col- 
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lapses to |+)b with probability cos^ or momentum 
—^Bi{Sn)\o) and spin state collapses to |—)_b with proba¬ 
bility sin^ When we obtain the outcome 'yBi{Sn)\o) 
there is no way to know if the spin state was rightly pro¬ 
tected or wrongly protected, as it is common to both 
cases. However when we obtain —jBi{Sn)\o), we come 
to know that the spin state was wrongly protected. In 
either of the cases we obtain no information about the 
unknown spin state, |m). Hence it’s impossible to clone 
by this method. 

Cloning in Center of Wave Packet(CWP) 
Frame: If we can somehow track the CWP of the split 
wave packets and measure its momentum, then we can 
perfectly clone an arbitrary unknown spin state, even 
with von-Neumann impulsive measurement, as follows: 
Consider the unknown state \rh) given in 8. Switch on in¬ 
homogeneous magnetic field Bizk for an interval of time 
T. Spin and sptial states evolve under the interaction 
Hamiltonian, 

Hint = —fl.Bizk = —^S.Bizk = —jSmm.Bizk 

= -yHj cos 9mSmZ (64) 

where, Sm = S.m, and S'm|m) = ||m). We are working 
in the eigenbasis of the operator, Sm- State at time t = t 
is, 

2 

\i’{T,Pz)) = exp{-^BiCosOmSmZT)\'m)xiPz) 
a 

= \m)x{Pz - xBi cos 9m^T) (65) 

Hence CWP has gained z-component of momentum equal 
to jBi cos9m^T. If we can some how measure this mo¬ 
mentum, we come to know 9m and there is no collapse. 
By following a procedure similar to that described from 
eq.s 49 to 55, we can find out (f>m- However if we work 
in the eigenbasis of the operator S^, we can observe the 
splitting. Hint can be written in terms of Sz as, 

Hint = -p-Bizk = -"fS.B^zk = -^BiSzZ ( 66 ) 

and the state at time t = t is exactly as given in eq.lO, 
where we can observe splitting and hence entanglement. 


CONCLUSION 

In principle it is possible to do exact calculations with¬ 
out the approximations that we have made. Hence, we 
conclude that, in principle it is possible to clone a single 
arbitrary unknown quantum state of a spin-1/2 particle 
(an electron) with arbitrary precision and with proba¬ 


bility tending to one. Of course, in practice we cannot 
make dimensions of double well potential, a and 6, and 
the time interval T, arbitrarily large, even though possi¬ 
ble in principle. Hence, it reduces the precision and suc¬ 
cess probability of cloning. However we note that, even a 
few millimeters and seconds are like infinity on the atomic 
scale. For eg., in NMR we talk of adiabatic processes even 
though life time of spins that we prepare is only a few 
seconds [29]. Hence, to get a quantitative picture of pre¬ 
cision and success probability attainable in practice, one 
has to do thorough calculation considering all practical 
limitations. Our protocol can be easily generalized to ar¬ 
bitrary charged, spin-1/2 particle. If it is not electrically 
charged, then one can exploit its gravitational potential 
energy, in principle, as splitting of wave packet is inde¬ 
pendent of charge, but depends only on spin. Protocol 
can also be generalized to higher spin (> 1/2) particles. 
Experimental realization of our protocol is with in the 
reach of current technology. 

Non-orthogonal state discrimination is more feasible 
than cloning, as it requires only the knowledge of polar 
angle, 9m- Also, to discriminate between, say, |0) and 
!+)(= (|0) + |l))/v^)i it is not necessary to find 9m with 
arbitrary precision. So much error is allowed such that 
we can just unambiguously discriminate between |0) and 
|-|-). Here we note that, inspite of all practical constraints 
on a, b etc., if we can still unambiguously discriminate 
between |0) and |-|-) in finite time T (such that T is less 
than time interval corresponding to two space-like sepa¬ 
rated events) with success probability greater than, say, 
0.5, then our protocol may open up the doors to super¬ 
luminal communication. To check this possibility, one 
has to do thorough analysis considering all practical con¬ 
straints. We also note that, the limit T —t oo corresponds 
to discriminating between |0) and |-|-) with success proba¬ 
bility tending to one, and it is not possible to do superlu- 
mial communication in this limit, as T is greater than any 
finite time interval corresponding to two space-like sep¬ 
arated events. But to communicate superluminally, we 
donot require success probability tending to one. Hence, 
there seems to be a possibility still. 
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o-i</>rr 


cot ^(cos — 1) 


ic/)„ 
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